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Abstract To explore and enrich the molecular mechanisms of thyroid hormone receptors (TRs) in the metamorphosis 
of amphibians, the cDNA sequences of 7Ra and TR in Microhyla fissipes were cloned and characterized. TRa was 1 
706 bp in length with an open reading frame (ORF) of 1 257 bp encoding a predicted protein of 418 amino acids and 
TRB was 1 422 bp with an ORF of 1 122 bp encoding a predicted protein of 373 amino acids. Their protein sequences 
contained 4 conserved domains of the nuclear receptor superfamily with two highly conserved cysteine-rich zinc fingers 
in the DNA-binding domain, whereas TRB was 42 amino acids shorter in its A/B domain than TRa. Highly-conserved 
sequences and structures indicated their conserved functions during metamorphosis. 7Ra expression reached peak at 
12 h and then decreased from 12 h to 48 h. While dramatically up-regulated TRP was observed after exposure of T3 
within 24 h, and it was down-regulated from 24 h to 48 h. The expression pattern of TR is similar to that in the natural 
metamorphosis. Furthermore, tadpoles treated 24 h also resembled the climax of metamorphosis tadpoles and TRS 
expression had higher responsiveness than TRa to T3 in M. fissipes. These results suggest M. fissipes may serve as the 
model to assay environmental compounds on TH signaling disruption. 
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1. Introduction 


Amphibian metamorphosis is a complex process, in 
which various organs and tissues undergo dramatically 
remodeling to transform from larva to juvenile 
(Atkinson et al., 1998; Brown and Cai, 2007; Du et al., 
2017). Although complex, this development process 
is completely initiated and orchestrated by only one 
hormone, thyroid hormone (TH) (Buchholz et al., 2006). 
The TH signaling pathway in amphibians has been 
well studied: when TH is absent, the unliganded TR/9- 
cis-retinoic acid receptor (RXR) heterodimers recruit 
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corepressors to repress the transcription of downstream 
target genes; once liganded with TH, the TR/RXR 
heterodimers undergo conformational change and thus 
allow to recruit coactivators to activate the same group of 
downstream genes (Morvan-Dubois et al., 2008; Grimaldi 
et al., 2013; Zhao et al., 2016). Therefore, TRs play 
important roles in the metamorphosis by acting as ligand- 
dependent transcriptional factors. 

TRs, members of a large superfamily of nuclear 
receptors (NR), possess a similar domain structure as 
that found in the other NRs: an N-terminal A/B domain 
with binding sites for transcriptional coregulators, a 
central DNA-binding domain C (DBD) containing two 
“zine fingers” for target gene recognition, a D domain 
(hinge region) containing the nuclear localization 
signal, and a C-terminal ligand-binding domain E/ 
F (LBD) where thyroid hormone binds and activates 
the receptor (Chen et al., 2014). There are two closely 
related families of TRs called TRa and TRf in vertebrate 
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(Yaoita et al., 1990; Chen et al., 2014). TRa and TR# are 
differentially expressed in various tissues of different 
species (Kawahara et al., 1991). Particularly, the TRa 
mRNA increases throughout the premetamorphosis 
stage of tadpole development, and falls after the climax 
of metamorphosis to a lower level in frogs (Yaoita and 
Brown, 1990). The TR mRNA is barely detectable 
during premetamorphosis. In synchrony with the onset 
of endogenous TH synthesis by the thyroid gland, the 
level of TR mRNA rises in parallel with endogenous 
TH, reaching a peak at the climax of metamorphosis and 
drops after metamorphosis (Choi et al., 2015). Although 
TH signaling pathway has been well studied, functional 
of TRa and TRP during metamorphosis have not been 
clearly characterized. TRa-deficient tadpoles developed 
faster with smaller body size than their wild-type siblings 
suggesting that TRa played important roles in controlling 
the timing of Xenopus tropicalis metamorphosis (Choi 
et al., 2015; Wen and Shi, 2015). Furthermore, disrupted 
TRa had different effect on the development of larval 
and juveniles and the metamorphosis of different 
organs (Choi et al., 2017, Wen et al., 2017). Different 
from TRa-knockout tadpoles, significantly delayed tail 
regression, the reduction in olfactory nerve length and 
head narrowing by gill absorption were detected in TRJ- 
knockout tadpoles (Nakajima et al., 2017). The different 
relative abundance levels of TRa and TRf transcripts 
induced by T3 where the general pattern was TRa = TRB 
in R. catesbeiana, while TRa < TRB in Xenopus laevis 
(Veldhoen et al., 2014). TRB was highly expressed during 
metamorphosis in M. fissipes and X. laevis, but TRa 
showed especially low expression in M. fissipes, implying 
that TR£ is essential for initiating metamorphosis, at least 
in M. fissipes (Zhao et al., 2016). 

Microhyla fissipes is a small-sized anuran from the 
family Microhylidae suborder Neobatrachia (Figure 
1). Due to the special expression pattern of TRs in M. 
fissipes mentioned above, it is important to clone TRs 
and understand the molecular mechanism of them in 
regulating metamorphosis. Furthermore, because of 
its characteristics (including wide distribution, fast 
development, development in vitro, strong survivability, 
biphasic life cycle, small body size, diploid and 
transparent tadpoles) and being induced to metamorphose 
by exogenous TH, M. fissipes may be an ideal model to 
evaluate the possible effects of environmental compounds 
on the thyroid system (Liu et al., 2016). 

Therefore, it is necessary to characterize TRs 
and evaluate responsiveness to TRs agonist (such 
as 3,3',5-Triiodo-L-thyronine, T3) in M. fissipes. 


Figure 1 Photograph of Microhyla fissipes. 


In this study, we have isolated, characterized, and 
phylogenetically analyzed TRa and TRf gene in 
M. fissipes, examined their expression pattern after T3 
treatment and explored the utility of M. fissipes as a model 
species for assaying TH signaling disrupting effects. 


2. Materials and Methods 


2.1. Animals sampling and experimental treatments 
Mature female and male M. fissipes were collected from 
Shuangliu, Chengdu, China (30.5825° N, 103.8438° 
E) in June, 2016. The male and female were injected 
luteinizing hormone-releasing hormone a (LHRHa) with 
0.3 ug/g body weight resolving dosage. Fertilized eggs 
were obtained from one pair of frogs and incubated in the 
dechlorinated tap water. Five days later, tadpoles were 
fed with spirulina powder once daily and subjected to a 
12:12 h light:dark cycle at 25 + 0.6°C. The developmental 
stage of tadpole was recorded using the M. fissipes 
developmental table (Wang et al., 2017). Tadpoles at 
stage 40 (metamorphosis climax) were selected for gene 


cloning. Tadpoles of stage 33 (premetamorphosis, oar- 
shaped limb bud) were selected to treat with 10 nmol/L 
3,3',5-triiodo-I-thyronine sodium (T3, Sigma-Aldrich, 
USA) for 48 h. The chemicals were renewed after 12 h of 
exposure when the medium was also refreshed. Tadpoles 
treated for 0 h, 12 h, 24 h, 36 h, and 48 h were collected 
for quantitative real-time (RT) PCR (n = 3 for each time 
point). Tadpoles treated 0 h were set as the control group. 
After anesthetization by MS222, tadpole sample was 
frozen immediately in liquid nitrogen, and then stored at 
—80°C for RNA extraction. 

The care and treatment of animals in this study were 
performed according to the Guideline for the Care 
and Use of Laboratory Animals in China. The animal 
experiments were approved by the Experimental Animal 
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Use Ethics Committee of the Chengdu Institute of 
Biology (Permit Number: 2016036). 


2.2. Cloning and molecular characterization of TRs 
Total RNA was extracted using TransZol (Transgen 
Biotech, Beijing, China), following the manufacturer’s 
instructions. Total RNA concentration was calculated 
using Nanodrop ND-1000 (Nanodrop, DW, USA). 
Partial cDNA sequence of TRs (TRa: comp77374_c0 ; 
TRB: comp124487_c3) were obtained from M. fissipes 
transcriptome (Zhao et al., 2016). According to the partial 
cDNA sequence, two specific primers for each gene, 
GSP3-1 and GSP3-2 (Table 1) were designed to amplify 
the 3' terminal regions by nested PCR. The 3' DNA 
ends were obtained using the SMART RACE cDNA 
Amplification Kit (Clontech, CA, USA) in accordance 
with the manufacturer’s instructions. Products of rapid 
amplification of cDNA ends (RACE) were cloned into 
pMD 18-T vector (TaKaRa, Japan) and sequenced 
using an automated DNA sequencer ABI3730 (Thermo 
Fisher Scientific, CA, USA) by Sangon Biotech Co. Ltd. 
(Shanghai, China). 

The amino acid sequence was deduced from the coding 
region via DNAStar (version 6.13). The cDNA sequence 
and the deduced amino acid sequence were analyzed 
using BLASTN and BLASTP, respectively. Deduced 
amino acid sequences of amphibian were aligned for 
analysis of putative conserved functional residues by 
Clustal X. The relevant amino acid sequences were 
obtained from the NCBI GenBank database: Rugosa 
rugose TRa BAM15695.1, Pelophylax nigromaculatus 
TRa AGT55994.1, Rana chensinensis TRa AIA98429.1, 
X. laevis TRaA NP_001081595.1, X. laevis TRaB 
BAL70322.1, X. tropicalis TRa NP_001039261.1, 
X. tropicalis TRB NP_001039270.1, X. laevis TRPA 
NP_001090182.1, X. laevis TRBB NP_001081250.1, P. 
nigromaculatus TRB AGT55995.1, R. chensinensis TRP 
AIA98430.1. 

Multiple cDNA sequences of TRa and TRP from 17 
species represented vertebrate (Mammalia, Aves, Reptilia, 
Amphibia, Pisces) and invertebrate were used in the 
sequence alignment by Clustal X. A phylogenetic tree 
was constructed by using the maximum likelihood (ML) 
method with the MEGA 6 (Tamura et al., 2013), and the 
reliability of the tree was assessed by the bootstrap method 
with 1,000 replications. The gene accession numbers 
are: Alligator mississippiensis TRa NM_001287278.1, 
Branchiostoma lanceolatum TRa EF672345.1, 
Crassostrea gigas TRa KP271450.1, Gallus gallus 
TRa NM_205313.1, Homo sapiens TRa AB307686.1, 
Oryzias latipes TRa AB114860.1, Rattus norvegicus 
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TRa M18028.1, P. nigromaculatus TRa KC139354.1, 
R. catesbeiana TRa L06064.1, R. chensinensis TRa 
KJ579109.1, R. rugose TRa AB683466.1, Schistosoma 
mansoni TRa AY395038.1, S. mansoni TRB AY395039.1, 
X. laevis TRaA M35343.1, X. laevis TRaB AB669465.1, 
X. tropicalis TRa NM_001045796.1, Nanorana 
parkeri TRa XM_018569566.1, Oryzias latipes TRa 
NM_001104705.1, Danio rerio TRa NM_131396.1, 
R. rugosa TRB AB683467.1, R. chensinensis TRB 
KJ579110.1, R. catesbeiana TRP L27344.1, P. 
nigromaculatus TRB KC139355.1, N. parkeri TRB 
XM_018570223.1, X. laevis TRBA NM_001096713.1, 
X. laevis TRGB NM_001087781.1, X. tropicalis 
TRP AB244214.1, G. gallus TRB NM_205447.2, A. 
mississippiensis TR NM_001287311.2, H. sapiens TRB 
M26747.1, R. norvegicus TRP J03933.1, O. latipes TRB 
AB114861.1, D. rerio TRB NM_131340.1. 


2.3. RNA isolation and gene expression analysis For 
TRs mRNA expression analysis after exposure of T3, 
we conducted qRT-PCR. Total RNA was extracted from 
tadpoles with TranZol reagent and first strand cDNA was 
synthesized from the same amount of RNA (1 ug) for 
each sample via TranScript All-in-One First-Strand cDNA 
Synthesis SuperMix for qPCR (TransGen, Beijing, China) 
with oligo (dT) primer. The expression of 7Ra and TRG 
mRNA was analyzed by qRT-PCR with the primer TRa 
and TRP (Table 1). Rp137 gene was used as a reference 
gene to normalize mRNA expression of 7Rs. PCR was 
performed in a reaction volume of 20 ul containing 10 ul 
of TransStart® Tip Green qPCR SuperMix (2x), 0.4 ul 
of Passive Reference Dye (50x), 0.4 ul each of forward 
and reverse primer (10 pmol/L), 8.2 ul of ddH,O, and 1 
ul of cDNA. Amplification was carried out in 7300plus 
(ABI, CA, USA), including 5 min at 95°C and 45 cycles 
of 5 s at 95°C and 31 s at 60°C, followed by a melting 
curve analysis. Each sample was run in triplicates. Each 
reaction was verified to contain a single product of the 
correct size by agarose gel electrophoresis. Quantitative 


data were shown as mean + SD (n = 3). The fold change 
of TRs expression after T3 treatment was determined 
by 2“ (cycle threshold, Ct). Data were then subjected 
to one-way analysis of variance (ANOVA) with SPSS 
Statistics 13.0 (SPSS Inc., Chicago, IL, USA). P < 0.05 
was regarded as statistically significant. 


3. Results 


3.1. Molecular characterization of TRs in M. fissipes 
The full-length cDNA sequences of TRa and TRJ were 
obtained by RNA-seq and 3'-RACE strategies. The full- 
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Table 1 Primers used in this study. F and R denote forward and reverse primer, respectively. 


Funticon Primers Sequence (5' to 3') Size (bp) Tm (°C) 
Cloning TRaGSP3-1 CCATCGCAAACACAACATTCCCCA 621 50 
TRaGSP3-2 TGACTTGCGTATGATCGGAGCCTG 567 60 
TRBGSP3-1 GCTGATGAAAGTCACCGACCTCCG 424 50 
TRBGSP3-2 CGCCAGCAGGTTCTTGCACATG 386 60 
3' CDS primer A AAGCAGTGGTATCAACGCAGACTAC(T) 
30 VN(N =A, C, G, or T; V =A, G, or C) 
CTAATACGACTCACTATAGGGCAAGCAGTGGTATCAACGCAGAGT 

a CTAATACGACTCACTATAGGGC 
RT-PCR TRa F: GCAGCCGTGCGCTATGAT 57 60 

R: T@CCATCTCACCGCTTAGTG 
TRB F: TGGCCAAAACTGCTGATGAA 56 60 

R: GCGTGGCAGGCTCCAA 

Rpl37 F: CCAAAAAGCGCAACAACCA 59 60 


length of TRa cDNA was | 706 bp in length and contained 
an open reading frame (ORF) of 1 257 bp, which encoded 
a peptide of 418 amino acids (Figure 2a). The TR cDNA 
was | 422 bp with an ORF of 1 122 bp, which encoded 
a peptide of 373 amino acids (Figure 2b). And the 3' 
untranslated region of TRa and TRB were 449 bp and 300 
bp, respectively. Two sequences were submitted to the 
GenBank (GenBank accession numbers: MG596879 and 
MGS596880). The homologies of nucleotide sequences 
and deduced amino acid sequences between 7Ra and 
TRP in M. fissipes were 61% and 72%, respectively. The 
calculated molecular weight of TRa polypeptide was 47.7 
kDa, and the theoretical isoelectric point (pI) was 7.08, 
while the calculated molecular weight of TRB polypeptide 
was 42.4 kDa with pI 6.76. 

A multiple alignment of TRs deduced amino acid 
sequences was performed in amphibians (Figure 3). 
The deduced amino acids of both TRa and TRB were 
composed of the N-terminal hypervariable region (A/ 
B domain), a conserved DNA-binding domain (DBD 
domain), a hinge region (D domain), and a ligand-binding 
domain (LBD domain). TRa and TRB in M. fissipes had 
high similarity in the DBD, D and LBD domains, whereas 
there was a deletion of 42 amino acids in TRB compared 
with TRa in the A/B domain (Figure 3). These deletion 
between TRa and TRB in M. fissipes corresponded 
to the differences found in X. laevis, X. tropicalis, R. 
chensinensis and P. nigromaculatus. Although there 
were several different sites in the A/B domain among 
TRs, the N-terminal signature sequence (NTSS) was 
well conserved. The conserved cysteine residues that 
comprise the two zinc fingers and the regulatory elements 
P-box as well as the T-box and A-box in the DBD were 
conserved. Furthermore, the consensus motif I (spanning 
helix 3—6) and motif II (from helix 8 to helix 10) and 


R: TTGCGAATCTGACGGACTTG 


the putative AF2 activation domain core (helix 12) were 
also identified. Therefore, this highly-conserved feature 
is likely to indicate pivotal significance of TRs in thyroid 
signaling pathway. 


3.2. Evolutionary relationships of TRs BLASTP and 
BLASTN in NCBI indicated that TRa and TRf shared 
different levels of homology with other species. The 
amino acid sequences of M. fissipes TRa had highest 
homolgy with R. rugosa, N. parkeri, R. chensinensis 
(98%, 10 different sites out of 418 amino acids), 
while M. fissipes TRB had highest homolgy with R. 
catesbeiana (98%, 8 different sites out of 373 amino 
acids). Nucleotide sequences of M. fissipes TRa and 
TRP shared high homology with N. parkeri (92% and 
93%), P. nigromaculatus (92% and 92%), R. rugosa 
(93% and 92%), X. tropicalis (88% and 87%), X. laevis 
A (88% and 80%), X. laevis B (88% and 78%) and had 
a lower homology with G. gallus (65% and 74%), A. 
mississippiensis (63% and 63%), H. sapiens (58% and 
47%), D. rerio (57% and 49%), B. lanceolatum (47% and 
49%), respectively. 

A phylogenetic tree constructed by the ML method 
from a multiple alignment of nucleotide sequences of M. 
fissipes TRa and TRP and a wide range of counterparts 
in various species including invertebrate, reptiles, 
birds, mammals and other amphibian (Figure 4). The 
phylogenetic tree showed that TRa and TR grouped into 
two highly consistent and separate clades. In both TRa and 
TRB clades, TRs from M. fissipes have similar positions 
in the phylogenetic tree. Furthermore, phylogenetic tree 
constructed based on TRs cDNA sequences was consistent 
well with the taxonomic positions of these organisms. 


3.3. Responsiveness of M. fissipes TRs expression 
to exogenous T3 In the presence of exogenous T3, 
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ATGGACCAGAATCTTAGCGGGCTGGACTGCTTGTCAGAGCCAGATGAAAAAAGGTGG 57 
M D Q NL S G LD C LSE P DE K R W 
CCGGATGGGAAGCGTAAAAGAAAGAACAGCCAATGCATGGGAAAAAGCGGCATGTCA 114 


P D G K RK R K N SRCMGKS GMS 
GGTGACCGTTCGGTGTCTCTGCTCTCTGCAGGGTACATCCCTAGCTACCTGGACAAAGAT 174 


G DRS V S L L S A G Y IPS YLDKD 
GAGCCGTGCGTGGTGTGCAGTGACAAGGCCACAGGGTACCACTACCGATGTATTACC 231 
E PC VV CS DK AT GY H YY R CIT 
TGTGAGGGTTGCAAGGGCTTCTTCCGTCGAACCATCCAGAAGAACCTGCACCCTTCC 288 


C E GC K G F F R R TIT QK N LH P S 
TACTCCTGCAAGTATGATGGCTGCTGCATCATCGATAAGATCACCCGGAATCAGTGCCAG 348 
Y S C K Y DG cCctteEdDK I TR N Q C Q 
CTCTGCCGATTCAAGAAATGTATCGCAGTGGGCATGGCAATGGACCTTGTCCTGGATGAT 408 
L C R F K K C I A V GM AM OD LV LDD 
TCCAAGCGGGTAGCCAAGCGGAAGCTGATTGAGGAAAATCGGGAACGGCGTCGGAAG 465 
S K R V A K R K LI E EN R E R RR K 
GAGGAGATGATGAAGTCTCTGCAGCAGCGACCTGAGCCAAGCAGCGAGGAATGGGAG 522 
E EM MK S L Q Q R P E PS SEE WE 
CTGATCCGCATAGTGACAGAAGCTCACCGGAGCACCAACGCACAGGGCAGCCACTGG 579 
L I RI V T EA H R S T N A Q GS H W 
AAGCAGCGCAGGAAGTTCTTGCCGGATGACATTGGGCAGTCTCCGATGGCCTCCATGCCA 639 
K Q RR K F L POD DI GQ S&S PM A S M P 
GATGGTGATAAAGTAGACCTGGAGGCGTTCAGTGAGTTCACCAAAATCATCACCCCTGCA 699 
D G D K V DL E A F § E F T K I I T PA 
ATCACCCGAGTGGTGGACTTCGCCAAGAAGCTGCCCATGTTCTCTGAGCTGCCCTGTGAA 759 
I T R V VD FA K K LPM F § EL P CE 
GACCAGATCATCCTGTTGAAAGGATGCTGTATGGAAATTATGTCTCTTCGGGCAGCCGTG 819 
D Qt1 1 %LLEK GE CM EITM § LR A A V 
CGCTATGATCCGGACAGCGAAACTCTGACACTAAGCGGTGAGATGGCAGTGAAGCAA 876 
R Y D PD S E T LT LS GEM A V K Q 
GAACAGCTGAAGAATGGAGGGCTGGGGGTTGTCTCTGATGCTATATTTGACTTGGGAAGA 936 
E QLKN GG LGVwVSs DAT FDL GR 
TCTTTGTCTGCGTTCAATCTGGATGACACAGAAGTGGCCCTGCTGCAGGCTGTCCTGCTC 996 
S LS A F N LDOD TT EV AL LQAWV LL 
ATGTCCTCAGACCGGACGGGTTTGATCTGCACGGACAAAATAGAGAAATGTCAGGAG 1053 


MS S DR TGLtI CC T DK I EK C QE 
ACGTACCTCCTCGCCTTCGAACACTACATCAACCATCGCAAACACAACATTCCCCACTTC 1113 
T ¥ EL &b AF BE H ¥ IN HR KE ANI PSE 
TGGCCCAAGCTCCTTATGAAAGTGACTGACTTGCGTATGATCGGAGCCTGCCATGCCAGT 1173 
w PK LLM K V T DLRM I GACHAS 


CGCTTCCTGCACATGAAGGTTGAGTGCCCCACAGAGCTCTTCCCGCCTCTCTTCCTTGAG 1233 
RFLHMKVECPTELFPPLFELE 
GTGTTTGAGGACCAGGAAGTCTGAGGGGCAGTGGAGTGGGCCAGAGACTAAAGGTAAT 1291 
V FE DQEV~- 
TITGCTTGGCTCTTTCCAGCCTGCCAACCCTTACAAGGCTCTTGIGICCTGCGAATGATG 1349 
CAAGGGAGTGGCGGGATAGGTCCTTCTGATAATGAAAACTGCCTTACGTAGTTGCAATG 1408 
CCTGGCCTTCCTGCTGCCTGATTTCTAACCAAGCTGGAGAGAGTTIGAAGCTITGCCACAG 1467 
CGCCTATCCTTCAGCAAAACTGTCCTTTITGACAGCATTAGTAAAGCACAGTGGTACTTA 1527 
AAGCAATAGATTCCATCTGGGCTCTTGTCGAGTCACTGTATGTAACATGGCCTGCTGAGA 1587 
GGAGCATGTCCTTAAAGCAGACCTTAACTCCCAGGGTAAACTGGCTATGTAACCCCCCA 1646 
CCAAGGACTTGAGTTAAGGGAAGTGGTGAAAAAAAAAAAAAAAAAAAAAAAAAAAA 1702 
AAAA 1706 


17 
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Continued Figure 2 


b. 

ATGCCTAGCAGCATGTCAGGGTACATCCCGAGCTACTTAGATAAAGATGAGCTGTGTGTG 60 
M P S S MS G Y I P S Y LD K D E L CV 
GTGTGTGGGGACAAGGCTACAGGATATCACTATAGATGCATCACCTGTGAGGGCTGCAAG 120 
V C G D KAT G Y HY R C I T C E GCK 
GGCTTTTTCCGGAGAACTATACAGAAAAATCTTCACCCTAGTTATTCATGTAAATATGAA 180 
G FF RR T I Q K N LHP S Y S C K YE 
GGAAAATGCGTGATAGACAAAGTAACCAGGAACCAGTGCCAAGAATGTCGCTTCAAA 237 
G KCV I DK V T R N Q CC Q E C RF K 
AAGTGTATTGCGGTTGGCATGGCAACAGACCTGGTTTTGGATGACAGCAAGCGTCTG 294 
K C I A V GM A T DL VL DOD S KRL 
GCAAAAAGAAGACTTATAGAAGAAAATCGAGAAAAGAGGCGCAAAGATGAAATACAG 351 
A KR R L I E ENRE KR RKDEI Q 
AAGTCGATGGTACAAAAACCCGAACCAACCCCGGAGGAGTGGGAGCTGATCCAGATT 408 
K SMV QKPEPTPEEWE L IQI 
GTCACCGAAGCTCACGTGACCACCAATGCACAGGGAAGTCACTGGAAGCAGAAACGG 465 
V T E A H V T T N AQ GS H WK Q K R 
AAATTTCTGCCAGAAGATATTGGGCAGGCTCCTATAGTTAATGCACCGGAGGGCGGAAAA 525 
K FL PEDI GQ A PT VNAPE GGK 
GTGGATCTAGAAGCTTTCAGCCAGTTTACAAAAATTATCACCCCAGCAATTACAAGAGTG 585 
V D L EA F § Q FT K I IT PA ITR V 
GTCGATTTTGCCAAAAAATTACCTATGTTTTGTGAGCTGCCATGTGAAGATCAGATCATC 645 
V DF AKKLPMFCE LP C E D QT I 
CTTCTTAAAGGCTGTTGCATGGAGATCATGTCGCTGCGAGCAGCAGTGCGCTACGACTCG 705 
L LK G CC ME I MS LRA AV RYDS 
GAAAGTGAGACCTTAACACTGAATGGGGAGATGGCAGTGACGCGAGGACAGCTAAAA 762 
E SE TLTLNGE MA V TRGOQLK 
AATGGAGGACTAGGAGTGGTATCTGACGCCATATTCGATTTGGGTGTTTCGCTGTCATCA 822 
N G G L G V V S D A I F DLG V S LSS 
TTTAATCTTGACGATACGGAAGTCGCCTTGTTGCAATCTGTGCTCCTAATGTCATCAGAT 882 
F NLDDT E VA LE Qs V LLMSSD 


CGTCCCGGTCTCTCCAGTGTGGAAAGAATAGAAAAGTGCCAAGAAGGTTTCCTCTTG 939 
R P GL SS V E RI E K CC Q E G FL L 
GCTTTTGAACACTACATTAACTACAGGAAACACAACGTTGCACACTTTTGGCCAAAA 996 
A FE H Y IN YR KH N V A H F W P K 
CTGCTGATGAAAGTCACCGACCTCCGCATGATTGGAGCCTGCCACGCCAGCAGGTTC 1053 
L LM KV T D LRM IGAC HAS REFE 
TTGCACATGAAAGTGGAGTGCCCCACTGAACTGTTCCCCCCACTGTTCTTAGAAGTG 1110 
L HM K V E C PTELFPP L FL E V 
TTTGAGGATTAGAAGAGACTGAGCTTCAGGATTCTCGGCACTACTGGGGTCTCCCCT 1167 
F E D - 


TCCATTCCATTGCCTCCGTCACTTTTTTGTTCTCGTCCGTTCTGAATAGACATAGATGA 1226 
AAATGTTCCTCTGATGCGGGTACTCGTGACTATCGTATTTTTGTTATTCTGTCCTTTTGA 1286 
TGTGAAGTTTATCCGTGGAAGTTTATCCCCTTTTGTGCTGAAGTGGCTGTCGACACATT 1345 
ACTCATCAGAGGGTTAGGCAGTCCTTCCAGTTATAAAGTGATTATTGTTTAAAAAAAA 1403 
AAAAAAAAAAAAAAAAAAA 1422 


Figure 2 cDNA and deduced amino acid sequences of Microhyla fissipes thyroid hormone receptor gene. (a) TRa; (b) TRG. The start codon 
(ATG) and the stop codon (TAG or TGA) are in bold. 
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5 15 25 35 45 55 15 85 95 105 115 125 135 145 155 165 
Rugosa rugosa TRa MDQNLSGLDC LSEPDEKRWP DGKRKRKNSQ CMGKSGMSGD SSVSLLS-— -AGYIPSYLD KDEI ICSD KATGYHYRCI TCEGCKGFFR RTIQKNLHPS YS¢KYDGCCI RFKKCIAV ILVLDD SKRVAKRKLI EENRERRRKE EMIKTLQQRP 
Pelophylax nigromaculatus TRa  MDQNLSGLDC LSEPDEKRWP DGKRKRKNSQ CMGKSGMSGD SSVSLLS—— -AGYIPSYLN KDEPQVVCSD KATGYHYRCI TCEGCKGFFR RTIQKNLHPS YS¢KYDGCCI RFKKCIAV ILVLDD SKRVAKRKLI EENRERRRKE EMIKTLQQRP 
Rana chensinensis TR MDQNLSGLDC LSEPDEKRWP DGKRKRKNSQ CMGKSGMSGD SSVSLLS—— -AGYIPSYLD KDEPQVVCSD KATGYHYRCI TCEGCKGFFR RTIQKNLHPS YS¢KYDGCCI RFKKCIAV ILVLDD SKRVAKRKLI EENRERRRKE EMIKTLQQRP 
Xenopus laevis TRA MDQNLSGLDC LSEPDEKRWP DGKRKRKNSQ CMGKSGMSGD SLVSLPS-—- -AGYIPSHLD KDEP@VVCSD KATGYHYRCI TCEGCKGFFR RTIQKNLHPS YS(KYDGCCI RFKKCIAV GMAMDLVLDD SKRVAKRKLI EENRQRRRKE EMIKTLQQRP 
Xenopus laevis TRa B MDQNLSGLDC LSEPDEKRWP DGKRKRKNSQ CMGKSGMSGD SLVSLPP-—- -AGYIPSYLD KDEP@VVCSD KATGYHYRCI TCEGCKGFFR RTIQKNLHPS YS@KYDGCCI RFKKCIAV GMAMDLVLDD SKRVAKRKLI EENRVRRRKE EMIKTLQQRT 
Xenopus tropicalis TRa MDQNLSGLDC LSEPDEKRWP DGKRKRKNSQ CMGKSGMSGD SLVSLPS—— -AGYIPSYLD KDEP@VVCSD KATGYHYRCI TCEGCKGFFR RTIQKNLHPS YS¢KYDGCCI RFKKCIAV GMAMDLVLDD TKRVAKRKLI EENRERRRKE EMIKTLQQRP 
Microhyal fissipes TRO MDQNLSGLDC LSEPDEKRWP DGKRKRKNSQ CMGKSGMSGD RSVSLLS—— -AGYIPSYLD KDEPGVVCSD KATGYHYRCI TCEGCKGFFR RTIQKNLHPS YSK RFKKCIAV GMAMDLVLDD SKRVAKRKLI EENRERRRKE EMMKSLQQRP 
Xenopus laevis TRBA MPSSMSVRLF TASAAQR--- —KKIQEGDC CVVLAGKTQG RFILIGAVAR VSGYIPSYLD KDEI (OGD KATGYHYRCI TCEGCKGFFR RTIQKNLHPS YS¢KYEGKCV IDKVTRNQ@CQ ECRFKKCKTV GMATDLVLDD SKRLAKRKLI EENREKRRKD EIQKSIVQRP 
Xenopus laevis TRBB MPSSMSVRLF TASAAQR——- —KKIQEGDC CVVLAGKTQG RFILIGAVAR VSGYIPSYLD KDEI (CGD KATGYHYRCI TCEGCKGFFR RTIQKNLHPS YS¢KYEGKCV IDKVTRNQCQ ECRFKKCKTV GMATDLVLDD SKRLAKRKLI EENREKRRKD EIQKSIVQRP 


Microhyal fissipes TRB MPSSMS- ad (OGD KATGYHYRCI TCEGCKGFFR RTIQKNLHPS YS¢KYEGKCV IDKVTRN@CQ ECRFKKCIAV GMATDLVLDD SKRLAKRRLI EENREKRRKD EIQKSMVQKP 
Pelophylax nigromaculatus TRB — MPSSMS- (CGD KATGYHYRCI TCEGCKGFFR RTIQKNLHPS YS¢KYEGKCV IDKVTRNQ@CQ ECRFKKCIAV GMATDLVLDD SKRLAKRKLI EENREKRRKD ELRKSLVQKP 
Rana chensinensis TRB MPSSMS- (OGD KATGYHYRCI TCEGCKGFFR RTIQKNLHPS YS¢KYEGKCV IDKVTRNQCQ ECRFKKCIAV GMATDLVLDD SKRLAKRKLI EENREKRRKD ELQKTLVQKP 
Xenopus tropicalis TRB MPSSMSETEE TKPG- (CGD KATGYHYRCI TCEGCKGFFR RTIQKNLHPS YS¢KYEGKCV IDKVTRNQ@CQ ECRFKKCIAV GMATDLVLDD SKRLAKRKLI EENREKRRKD EIQKSLVQKP 
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Rugosa rugosa TRa EPSSEEWELI RIVTEAHRST NAQGSHWKQR RKFLPEDIGQ NPMASMPDSD KVDLEAFSÈF TKIITPAITR VVDFAKKLPM FSELPCEDQI ILLKGCCMEI MSLRAAVRYD HDSETLTLSG EMAVKREQLK NGGLGVVSDA |[FDLGRSLSA FNLDDTEVAL LQAVLLMSSD 
Pelophylax nigromaculatus TRa _ EPSSEEWELI RIVTEAHRST NAQGSHWKQR RKFLPEPIGQ NPMASMPDSD KVDLEAFSEF TKIITPAITR VVDFAKKLPM FSELPCEDQI ILLKGCCMEI MSLRAAVRYD ETLILSG EMAVKREQLK NGGLGVVSDA |IFDLGRSLSA FNLDDTEVAL LQAVLLMSSD 
Rana chensinensis TR EPSSEEWELI RIVIEAHRST NAQGSHWKQR RKFLPEDIGQ NPMASMPDND KVDLEAFSEF TKIITPAITR VVDFAKKLPM FSELPCEDQI ILLKGCCMEI MSLRAAVRYD HDSETLTLSG EMAVKREQLK NGGLGVVSDA [IFDLGRSLSA FNLDDTEVAL LQAVLLMSSD 
Xenopus laevis TRa A EPSSEEWELI RIVIEAHRST NAQGSHWKQR RKFLPEDIGQ SPMASMPDGD KVDLEAFSRF TKIITPAITR VVDFAKKLPM FSELTCEDQI ILLKGCCMEI MSLRAAVRYD HDSETLTLSG EMAVKREQLK NGGLGVVSDA |[FDLGRSLAA FNLDDTEVAL LQAVLLMSSD 
Xenopus laevis TRa B EPSSEEWELI RIVIEAHRST NAQGSHWKQR RKFLPEDIGQ SPMASMPDGD KVDLEAFSRF TKIITPAITR VVDFAKKLPM FSELTCEDQI ILLKGCCMEI MSLRAAVRYD ETLTLSG EMAVKREQLK NGGLGVVSDA |IFDLGRSLAA FNLDDTEVAL LQAVLLMSSD 
Xenopus tropicalis TRa EPSKEEWELI RIVIEAHRST NAQGSHWKQR RKFLPEDIGQ SPMASMPDGD KVDLEAFS§F TKIITPAITR VVDFAKKLPM FSELPCEDQI ILLKGCCMEI MSLRAAVRYD HDSETLTLSG EMAVKREQLK NGGLGAVSDA |IFDLGRSLAA FNLDDTEVAL LQAVLLMSSD 
Microhyal fissipes TRa EPSSEEWELI RIVIEAHRST NAQGSHWKQR RKFLPDDIGQ SPMASMPDGD KVDLEAFSRF TKIITPAITR VVDFAKKLPM FSELPCEDQI ILLKGCCMEI MSLRAAVRYD HDSETLTLSG EMAVKQEQLK NGGLGVVSDA |I[FDLGRSLSA FNLDDTEVAL LQAVLLMSSD 
Xenopus laevis TRBA EPTQEEWELI QVVTEAHVAT NAQGSHWKQK RKFLPEDIGQ APIVNAPEGG KVDLEAFSQF TKIITPAITR VVDFAKKLPM FCELPCEDQI ILLKGCCMEI MSLRAAVRYD HESETLTLNG EMAVTRGQLK NGGLGVVSDA |IFDLGVSLSS FSLDDTEVAL LQAVLLMSSD 
Xenopus laevis TRBB EPTQEEWELI QVVTEAHVAT NAQGSHWKQK RKFLPEDIGQ APIVNAPEGG KVDLEAFSQF TKIITPAITR VVDFAKKLPM FCELPCEDQI ILLKGCCMEI MSLRAAVRYD HESETLTLNG EMAVTRGQLK NGGLGVVSDA |IFDLGVSLSS FSLDDTEVAL LQAVLLMSSD 
Microhyal fissipes TRB EPTPEEWELI QIVTEAHVTT NAQGSHWKQK RKFLPEDIGQ APIVNAPEGG KVDLEAFS@F TKIITPAITR VVDFAKKLPM FCELPCEDQI ILLKGCCMEI MSLRAAVRYD GESETLTLNG EMAVTRGQLK NGGLGVVSDA |[FDLGVSLSS FNLDDTEVAL LQSVLLMSSD 
Pelophylax nigromaculatus TRB — EPTSEEWELI QVVTEAHVAT NAQGSHWKQK RKFLPEDIGQ APIVNAPEGG KVDLEAFSQF TKIITPAITR VVDFAKKLPM FCELPCEDHI ILLKGCCMEI MSLRAAVRYD HESETLTLNG EMAVTRGQLK NGGLGVVSDA |IFDLGVSLSS FNLDDTEVAL LQAVLLMSSD 
Rana chensinensis TRB EPTSEEWELI QVVTEAHVAT NAQGSHWKQK RKFLPEDIGQ APIVNAPEGG KVDLEAFSQF TKIITPAITR VVDFAKKLPM FCELPCEDQI ILLKGCCMEI MSLRAAVRYD HESETLTLNG EMAVTRGQLK NGGLGVVSDA ILGVSLSS FNLDDTEVAL LQAVLLMSSD 
Xenopus tropicalis TRB EPTQEEWELI QVVSEAHVAT NAQGSHWKQK RKFLPEDIGQ APIVSAPEGG KVDLEAFSQF_TKIITPAITR VVDFAKKLPM FCELPCEDQI ILLKGCCMEI MSLRAAVRYD HESETLTLNG EMAVTRGQLK NGGLGVVSDA |LFDLGVSLSS_FSLDDTEVAL LQAVLLMSSD 
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Rugosa rugosa TRa RIGLICIDKI EKCQETYLLA FE 
Pelophylax nigromaculatus TRa RTGLICTDKI EKCQETYLLA FE 


YINHRKH NIPHFWPKLL MKVTDLRMIG ACHASRFLHM KVECPTELFP PLFLEVFEDQ EV 
INHRKH NIPHFWPKLL MKVIDLRMIG ACHASRFLHM KVECPTELFP PLFLEVFEDR EV 
INHRKH NIPHFWPKLL MKVIDLRMIG ACHASRFLHM KVECPTELFP PLFLEVFEDQ EV 
'INHRKH NIPHFWPKLL MKVTDLRMIG ACHASRFLHM KVECPTELFP PLFLEVFEDQ EV 
INHRKH NIPHFWPKLL MKVIDLRMIG ACHASRFLHM KVECPTELFP PLFLEVFEDQ EV 
YINHRKH NIPHFWPKLL MKVIDLRMIG ACHASRFLHM KVECPTELFP PLFLEVFEDQ EV 
YINHRKH NIPHFWPKLL MKVIDLRMIG ACHASRFLHM KVECPTELFP PLFLEVFEDQ EV 


Rana chensinensis TRa RIGLICTDKI EKCQETYLLA FE 
Xenopus laevis TRa A RIGLICTDKI EKCQETYLLA FE 
Xenopus laevis TRaB RIGLICTDKI EKCQETYLLA FE 
Xenopus tropicalis TRa RIGLMCTDKI EKCQETYLLA FE 
Microhyal fissipes TRa RIGLICTDKI EKCQETYLLA FE 
Xenopus laevis TRBA RPGLSSVERI EKCQEGFLLA FEHYINYRKH NIAHFWPKLL MKVIDLRMIG ACHASRFLHM KVECPTELFP PLFLEVFED- — 
Xenopus laevis TRBB RPGLSSVERI EKCQEGFLLA FE 
Microhyal fissipes TRB RPGLSSVERI EKCQEGFLLA FE 
Pelophylax nigromaculatus TRB — RPGLSSVERI EKCQEGFLLA FE 


YINYRKH NIAHFWPKLL MKVTDLRMIG ACHASRFLHM KVECPTELFP PLFLEVFED- — 
YINYRKH NVAHFWPKLL MKVTDLRMIG ACHASRFLHM KVECPTELFP PLFLEVFED- — 
INYRKH NVAHFWPKLL MKVIDLRMIG ACHASRFLHM KVECPTELFP PLFLEVFED- — 
Rana chensinensis TRB RPGLSSVERI EKCQEGFLLA FEHYINYRKH NVAHFWPKLL MKVTDLRMIG ACHASRFLHM KVECPTELFP PLFLEVFED- — 


Xenopus tropicalis TRB RPGLTSVERI EKCQEGFLLA FEHYINYRKH NIAHFWPKLL MKVTDLRMIG ACHASRFLHM KVECPTELFP PLFLEVFED- — 
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Figure 3 Alignments of TR amino acid sequences from Microhyla fissipes with other amphibian species. Asterisk (*) indicated conserved 
amino acids and hyphens (—) represented spaces inserted to maximize similarity. Colon (:) indicates conservation between groups of strongly 
similar properties, and period (.) indicates conservation between groups of weakly similar properties. The black vertical lines indicate the 
borders of four domains: the N-terminal hypervariable region (A/B domain), DNA-binding domain (DBD domain), hinge region (D domain) 
and ligand binding domain (LBD domain). Triangle indicated the conserved cysteine residues that comprise the two zinc finger with red box 
of the DBD. The conserved N-terminal signature sequence (NTSS), DR-box, P-box, T-box, A-box and Helices (H3-H12) are figured out. 
Motif I and Motif II are boxed by blue while activation domain (AF2-AD) is purple double underlined. 


the morphology of M. fissipes tadpoles has changed increased by T3 at 12 h and 24 h, respectively. And then 
dramatically, and the characteristic of morphology change TRa and TRJ expression decreased to the lower level 
trend was similar to that during the natural metamorphosis at 48 h than 0 h (Figure 5). Furthermore, the same gene 
(data not shown). To determine whether exogenous T3 expression pattern of TR has been detected between 
could potentially regulate the expression pattern of TRa exogenous T3 induced and natural metamorphosis (Zhao 
and TRG, we exposed the premetamorphic tadpoles to T3 etal., 2016). Therefore, different responsiveness of TRa 
and analyzed by qPCR. These results showed that M. and TR£ to T3 indicated different functions of them in the 
fissipes TRa and TRP mRNA expression were significantly | metamorphosis. 
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Figure 4 Phylogenetic tree of nucleotide sequences of TRs genes of different species using the the maximum likelihood method. Numbers at 
the branches represent the bootstrap support values. Microhyla fissipes TRa and TRf are highlighted in the box, respectively. 


4. Discussion 


identities. The expression patterns after T3 treatment 
indicated the important roles of TRa and TR/ during the 


This study was designed to determine TRa and TR 
sequences and analyze their expression patterns after T3 
exposure to gain further insights as to how these genes 
may function in metamorphosis. Therefore, TRa and 
TRB genes of M. fissipes were cloned by RNA-seq and 
RACE. Phylogenetical analysis showed M. fissipes TRa 
and TRJ gene had high homology with the corresponding 
genes of other amphibians at both the nucleotide and 
the amino acid level, respectively, confirming their 


metamorphosis. 

Both TRa and TRB amino acid structure identified in 
M. fissipes possessed the typically functional domains 
of the NR superfamily. All of them had a conserved 
DBD domain with two zinc fingers, which determined 
the binding specificity of nuclear receptors (Natalia 
and Thorsten, 2004), and a LBD domain containing the 
typical 12 helices (Marchand et al., 2001). The P-box 
determining DNA binding specificity interacted with 
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Figure 5 Relative expression levels of TRs genes after exposure 
to 10 nmol/L T3. Data were expressed as the mean fold difference 
(mean + SD, n = 3) from the control group. Significant differences 
between treatment groups and the controls were indicated by 
* (P <0.05). 


the specific response element AGGTCA, and the T-box 
and A-box regions contributed to dimerization and 
DNA binding stabilization, respectively (Manchado et 
al., 2009). The conserved NTSS (GYIPS(Y/H) L(D/N) 
KDE(P/L)) which was the TR specific motif was also 
detected in the C-terminus of the variable A/B domain 
of M. fissipes TRs (Wu et al., 2007). The deletion of 
42 amino acids in the A/B domain of TRB indicated its 
different function from TRa. The conserved AF2-AD 
motif (LFLEVF) played an important role in recruiting a 
coactivator (Nagy et al., 1999; Nelson and Habibi, 2009). 
The conservation of structure and functional roles of the 
above mentioned sites in TRs would be consolidated by 
their high identity throughout the evolution of vertebrates. 

TRa and TRJ have been demonstrated to have high 
homology across vertebrates (Oka et al., 2013). In 
this study, the amino acid and nucleotide sequence 
homologies of 7Ra and TRJ between M. fissipes and 
other amphibians were over 90%, but lower homology 
with invertebrates. In the phylogenetic tree, vertebrate 
TRa and TRB were located in two clearly separated 
clades, in accordance with the fact that TRa and TR 
may be the products of an ancient gene duplication event 
during evolution (Chen et al., 2014). Furthermore, the 
homologies of TRa and TRf to the corresponding genes 
from other species accorded with their evolutionary 
relationship. Only one 7Ra and TRJ gene were identified 
in amphibian except in X. laevis which is tetraploid. The 
homologies of nucleotide sequences and deduced amino 
acid sequences between TRa and TRP in M. fissipes were 
61% and 72%, respectively. The homology of TRa to 
TRP in M. fissipes was lower than R. nigromaculata (72% 
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and 86%) and X. laevis (74% and 85%). Low sequence 
identification of TRa and TR# may indicate their different 
regulated function in M. fissipes metamorphosis, which is 
also implied by their different expression pattern during 
natural metamorphosis (Zhao et al., 2016). 
Metamorphosis is a critical developmental stage 
mediated by TH in amphibian (Wang et al., 2008). 
The function of TRs as transducers of TH responses 
has converted NRs in targets to clarify the molecular 
mechanisms that govern metamorphosis (Manchado et 
al., 2009). To understand the potential importance of 
these two receptors in M. fissipes metamorphosis and 
function of exogenous ligand for the receptor systems, we 
examined their mRNA expression levels by qPCR after T3 
exposure. T3 not only induced M. fissipes premetamorphic 
tadpoles metamorphosis at the morphology and histology 
level (data not shown), but also induced TRa and TRB 
expression, which have also been reported in X. laevis, 
X. troplis and R. catesbeiana (Shi, 1999; Wang et al., 
2008). During the natural metamorphosis, TR expression 
increased dramatically and correlated with the endogenous 
THs, while TRa expression slightly increased; and all 
of them decreased at the end of metamorphosis (Shi, 
1999; Zhao et al., 2016). In this study, TRa expression 
reached peak at 12 h and then decreased from 12 h to 48 
h. While dramatically up-regulated TRJ expression was 
observed after exposure of T3 within 24 h, and it was 
down-regulated from 24 h and with the lowest expression 
observed at 48 h of T3 treatment. These results suggested 
that TRJ expression pattern after T3 treatment is the 
same as that during natural metamorphosis, and the 
expression of TR# in tadpoles treated with T3 for 24h also 
resembled its expression in the tadpoles at the climax of 
metamorphosis (Zhao et al., 2016), which correlated with 
the morphological and histological results. Therefore, 
we can use T3 to simulate metamorphosis for further 
research on metamorphosis which will be high-efficiency 
and time saving. Moreover, M. fissipes could also 
serve as the model to assay environmental compounds 
on TH signaling disruption, while expression of TRs 
(in particular TRJ) has been also used as a molecular 
biomarker for assaying TH signaling disrupting actions in 
X. laevis and R. nigromaculata because of their response 
to TH (Opitz et al., 2006; Veldhoen et al., 2006; Lou et 
al., 2014). Due to the specific expression level of TRs 
during M. fissipes natural metamorphosis and dramatically 
up-regulated TRS mRNA expression after T3 treatment, 
we will use genome editing tools such as CRISPR/Cas9 
to illustrate the mechanism of TRa and TRP in M. fissipes. 
In conclusion, TRa and TR from M. fissipes were cloned 
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and characterized for the first time. Functional site and 
phylogenetic analysis indicated the conserved function 
of TRs from invertebrate to vertebrates. Premetamorphic 
tadpoles treated with T3 for 24 h resembled the climax of 
metamorphosis tadpoles during natural metamorphosis, 
and TR mRNA expression analysis could serve as a 
sensitive molecular testing approach to study effects 
of environmental compounds on the thyroid system in 
M. fissipes. 
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